Spectral Representations of Neutron-Star Equations of State 
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Methods are developed for constructing spectral representations of cold (barotropic) neutron-star 
equations of state. These representations are faithful in the sense that every physical equation 
of state has a representation of this type, and conversely every such representation satisfies the 
minimal thermodynamic stability criteria required of any physical equation of state. These spectral 
representations are also efficient, in the sense that only a few spectral coefficients are generally 
required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency 
is illustrated by constructing spectral fits to a large collection of "realistic" neutron-star equations 
of state. 



I. INTRODUCTION 

The gravitational field of a neutron star compresses the 
material in its core to densities that exceed those inside 
normal atomic nuclei. The resulting matter is a mixture 
of free baryons (neutrons and protons), leptons (electrons 
and muons), and likely also smaller fractions of hyperons, 
mesons, or perhaps even free quarks. The basic thermo- 
dynamic relationship between the total energy density e, 
and the pressure p of this material is called its equation 
of state, e = e(p), and is determined by the complicated 
micro-physical interactions between the various particle 
species present in the mixture 

In addition to its dependence on the pressure, the en- 
ergy density of a mixture also depends typically on the 
temperature and other quantities like the relative abun- 
dances of the different particle species. Thus there is no 
guarantee that a simple barotropic form, e = e(p), applies 
to any, let alone universally to all, neutron-star matter. 
Yet there is reason to expect that a universal barotropic 
form might be an excellent approximation. Neutron stars 
are born when the cores of massive stars (and perhaps 
white dwarfs) become unstable and undergo gravitational 
collapse. Compression heats the material as it collapses, 
to temperatures that exceed the binding energies of all 
atomic nuclei. Neutron-star matter always begins there- 
fore as a very hot plasma of free baryons and leptons, 
etc. This material is expected to evolve quickly to the 
lowest available energy state as it cools by neutrino and 
photon emission, and this fixes the relative abundances of 
the various particle species. The thermal energies of the 
particles fall rapidly below their Fermi levels, so the ther- 
mal contribution to the energy rapidly becomes negligi- 
ble. The matter therefore is expected to evolve on a very 
short time scale to a state that is well described by a tem- 
perature independent barotropic equation of state. This 
paper develops more efficient ways to represent equations 
of state of this type. 

The matter densities in the cores of neutron stars are 
well beyond the reach of current laboratory experiments. 
Heavy-ion scattering (including most recently those con- 
ducted at RHIC and LHC) provides a wealth of informa- 
tion about the interactions among the various particles 
expected to make up neutron-star matter. Unfortunately 



those experiments bear only indirectly on the properties 
of the equilibrium ground state, because the effective 
temperature of the nuclear matter in the experiments 
is quite high. Little insight is provided therefore into 
effects, like complicated many body interactions, that 
might play a role only in states with low temperature 
and high density and pressure. 

There has also been a significant effort over the past 
several decades to understand this material from a theo- 
retical perspective: hundreds of papers devoted to model- 
ing neutron-star matter have appeared in the literature. 
But the properties of this material are far outside the 
realm where the usual arsenal of theoretical tools were 
designed to work reliably, so it is not surprising that there 
is no consensus among theoreticians yet on the neutron- 
star equation of state. For example, the current models' 
predictions of the pressure at a given density still vary 
by about an order of magnitude [2| . 

Direct observations of neutron stars may be the most 
promising approach to understanding the properties of 
high-density nuclear matter. It is well known that the 
equation of state along with the gravitational field equa- 
tions determine the observable macroscopic properties of 
neutron stars [|| , and conversely that a complete knowl- 
edge of an appropriate set of macroscopic properties (e.g. 
masses and radii) determines the equation of state [J]. 
Studies of neutron-star models show that their macro- 
scopic properties, like their masses and radii, vary widely 
even within the current "realistic" class of equations of 
state 0, [1| . So it has long been recognized that accurate 
observations of the macroscopic properties of neutron 
stars will provide significant constraints on the nuclear 
equation of state. Unfortunately the needed observations 
are quite difficult to make. Masses of several dozen neu- 
tron stars have now been measured fairly accurately (see 
e.g. Ref. 0), and these observations have ruled out large 
classes of very "soft" equation of state models. Only a 
few radius measurements have been made however Q, 
and these are not reliable and accurate enough yet to 
make solid quantitative measurements of the equation of 
state itself possible. 

There is reason to hope that more abundant and ac- 
curate measurements of both neutron-star masses and 
radii will become available, however. When the first ac- 
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curate measurements are made, they are not likely to 
be numerous and accurate enough to determine the en- 
tire high density portion of the neutron-star equation 
of state. Various attempts have been made, therefore, 
to find representations of equations of state that make 
their essential features depend on just a few parameters. 
One approach is to use the parameters that characterize 
the nuclear interaction models as a way to parametrize 
the equations of state constructed from them. These 
might include a number of micro-physical parameters like 
the bulk nucleon incompressibility and symmetry energy 
parameters in models of the nucleon interaction poten- 
tial , or the coupling constants and mixing angles in ef- 
fective mean field theory descriptions . Comparing the 
masses and radii of neutron stars based on these model 
equations of state should fix the values of the unknown 
nuclear interaction model parameters. This approach 
would clearly be ideal if a reliable and accurate micro- 
physical model of neutron-star matter were known. Un- 
fortunately there is no consensus that any of the existing 
nuclear-matter models are good enough yet to describe 
neutron-star matter accurately and reliably. 

Another approach is to construct purely empirical fits 
rather than micro-physics based models of the equation 
of state. The first attempts to do this [3, EI3 approxi- 
mated the high density part of the equation of state as a 
simple polytrope, i.e., an equation of state in which the 
pressure is proportional to a power of the density. 1 These 
first simple fits were shown to reproduce the central pres- 
sures and densities of neutron-star models based on "re- 
alistic" equations of state with about 15% accuracy [J]. 
This type of approximation can be improved by dividing 
the relevant range of pressures into a number of intervals, 
Po < pi < ... < p m ax, and fitting a different polytrope to 
the equation of state in each interval. Any level of ac- 
curacy can then be achieved by using a sufficiently large 
number of intervals. A number of authors have proposed 
using piccewisc polytropes to approximate the high den- 
sity part of the neutron-star equation of state @, [13, [H| > 
and these approximations turn out to be quite efficient. 
Fits of the high density parts of "realistic" neutron-star 
equations of state have been shown to achieve accura- 
cies of a few percent for piecewise polytropes with only 
a small number of free parameters [2], [UJ, LL!| ■ The most 
extensive study to date uses fits with four free param- 
eters that give average errors of only a few percent for 
34 realistic equations of state @] . Other types of empiri- 
cal fits to neutron-star equations of state have also been 
reported in the literature fl2T4l5l |. These provide high 
accuracy approximations of particular realistic equations 



The term relativistic polytrope is most commonly used for equa- 
tions of state that satisfy p oc p r , where p is the conserved rest- 
mass density and T is the (constant) adiabatic index. It is also 
used (less commonly) for those satisfying p oc e 1 where e is the 
total energy density and 7 is a constant. The two definitions 
agree in the Newtonian limit. 



of state (generally using fifteen to twenty parameters to 
do this), and do not appear to have been intended as 
efficient ways to model large classes of equations of state. 

This paper continues the effort to construct effi- 
cient empirical representations of "realistic" neutron-star 
equations of state. New methods are described here for 
constructing parametric representations based on spec- 
tral fits. Spectral representations are generalizations of 
the Fourier series used to represent periodic functions. 
It is shown in Sec. [TT] that spectral representations can 
be constructed that are faithful, in the sense that every 
physical equation of state has such a representation and 
conversely that every such representation satisfies the ba- 
sic thermodynamic stability conditions required of any 
equation of state. It is also shown in Sec. IIIII that these 
spectral representations do a good job of representing 
the currently available "realistic" neutron-star equations 
of state. A suitably constructed two-parameter spectral 
representation is shown to be about as accurate as the 
most carefully studied four-parameter polytrope fits @]. 
For smooth equations of state, the errors in the spectral 
fits decrease exponentially as the number of parameters 
is increased, while piecewise polytrope fits generally de- 
crease only quadratically. These spectral fits make it 
possible therefore to provide very accurate representa- 
tions of the neutron-star equation of state using only a 
small number parameters. They should provide an im- 
portant new tool for extracting the high density equation 
of state from neutron-star observations. 



II. SPECTRAL REPRESENTATIONS OF THE 
EQUATION OF STATE 

Any equation of state, e = e(p), can be represented in 
a "spectral" expansion, as linear combinations of basis 
functions $fc(p): 

e(p) = J2^k(p)- (1) 

k 

Any complete set of functions, such as the Fourier basis 
functions or the Chebyshev polynomials, could be used 
as the in these expansions. Equations of state are 
determined in such representations by their spectral co- 
efficients, £fc. Truncated versions of these expansions, in 
which only a finite number of terms are kept, provide ap- 
proximate parametric representations of arbitrary equa- 
tions of state: e = e(p, e^). 

Physical equations of state must be non-negative, 
e (p) > 0, and monotonically increasing functions, 
de{p)/dp > 0, to insure thermodynamic stability. Since 
almost all functions fail to satisfy these conditions, it fol- 
lows that almost all choices of spectral coefficients, in 
an expansion such as Eq. (TTJ represent functions that can 
not be equations of state. Thus the spectral coefficients 
£fe obtained by fitting to a physical equations of state 
are likely to produce a representation that violates basic 
thermodynamic stability. So unfortunately, representing 



3 



an equation of state with a straightforward spectral ex- 
pansion is not particularly useful. 

Instead, faithful representations are needed: ones that 
ensure the positivity and monotonicity conditions for ev- 
ery choice of spectral coefficients. Methods of construct- 
ing faithful representations of the equation of state are 
presented in the following sections. Spectral represen- 
tations of the standard form of the equation of state, in 
which the energy density is expressed as a function of the 
pressure e = e(p), are presented in Sec. Ill Al For some ap- 
plications it is more convenient to express the equation 
of state in terms of the relativistic enthalpy e = e(h). 
Spectral representations of these enthalpy based forms 
are given in Sec. Ill Bl 



A. Pressure Based Forms 

An equation of state e(p) determines, and is deter- 
mined by (up to an integration constant), the adiabatic 
index r(p), defined by 



r(p) 



pdp 



P 



de 



(2) 



Given the equation of state e(p) is determined sim- 
ply by integrating the first-order ordinary differential 
equation, 



de(p) _ e(p) +p 



dp 



pT(p) 



(3) 



The adiabatic index must be positive T(p) > to ensure 
thermodynamic stability, but it need not be monotonic. 
Thus a larger class of functions represent possible phys- 
ical adiabatic indices, and this makes it easier to rep- 
resent equations of state through spectral expansions of 
r(p). In particular, every physical equation of state can 
be represented by the following spectral expansion of the 
adiabatic index T(p): 



T(p) = cxp 



^lk^k(p) 



(4) 



Conversely, every choice of (for which the series in this 
expansion converges) results in a positive adiabatic index, 
and thus an equation of state from Eq. ([3]), that satisfies 
the positivity and thermodynamic stability conditions. 
So this representation is faithful in the sense defined in 
Sec. EH 

The construction of an explicit spectral representation 
of the equation of state requires a choice for the basis 
functions $k{p)- To that end, it is useful to define a 
dimcnsionless logarithmic pressure variable: 



x = log(p/p ). 



(5) 



The constant po is a scale factor, chosen here to be the 
minimum value of the pressure, po < p, in the domain 



where the spectral expansions are to be used. The fol- 
lowing expansion of the adiabatic index T(x) is found to 
be useful and effective, 



T(x) = exp 




(6) 



One advantage of this simple power-law basis is that the 
lowest order spectral coefficients 7^ have fairly simple 
physical interpretations. For example the lowest order 
coefficient, 70, is determined by the adiabatic index eval- 
uated at the reference pressure: 70 = logr(po)- Simi- 
larly, the next coefficient, 71, determines the behavior of 
the adiabatic index, T(p) sw r(po)(p/po) 71 f° r pressures 
near p , i.e. for x = log(p/p ) -C 1. 

The power law basis, $k{%) = % , has the advantage 
of simplicity. While it might be advantageous to choose 
another basis like the Chebyshev polynomials for some 
purposes, these advantages can only be fully exploited 
by using a knowledge of the exact range, < x < 
and re-scaling x in the optimal way. The additional in- 
formation, like a; max for example, needed to do that will 
not be available a priori for the real neutron-star equa- 
tion of state. So here the simple power law basis is used, 
and fortunately this choice seems to work quite well. 

Given an adiabatic index, T(p), it is straightforward 
to determine the equation of state, e(p), by integrating 
the ordinary differential equation, Eq. ([3]). The solutions 
and hence the equation of state can be reduced to quadra- 
tures: 



1 



fJ-ip) Hip) Jpo t (p') 
where fi(p) is defined as 



Hip) = exp 



dp' 



p'T(p') 



(7) 



(8) 



and where eo = e(po) is the constant of integration needed 
to fix the solution. This eo is fixed in the fits performed 
in Sec. IHII bv matching to a low density equation of state 
at the pressure, po (i.e. at xq = 0), chosen to be a point 
somewhat below nuclear density. 

The quadratures indicated in Eq. ([7]) can not be done 
analytically for the expansion given in Eq. (j6|), so an ex- 
plicit analytic expression for the equation of state is not 
available in this case. However, the integrands in these 
quadratures are analytic functions that can be integrated 
numerically very accurately and efficiently Using Gaus- 
sian quadrature for example, double precision accuracy 
can be achieved using about 10 points for each integral. 
So there is very little practical difference between having 
an explicit analytic expression for the equation of state, 
and the expression in Eq. ([7J in terms of quadratures of 
explicit analytic functions. 

It might be advantageous in some situations to con- 
struct spectral expansions for the equation of state using 
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thermodynamic quantities other than T(p). For example, 
the adiabatic sound speed, v(p), defined by 



v 2 {p) 



de(p) 



dp 



(9) 



(where c is the speed of light) could be used to obtain 
the equation of state by integrating the simple ordinary 
differential equation, 



de(p) c 2 
dp v 2 (p) 



(10) 



The thermodynamic stability condition, < v 2 , could 
be enforced in this case by constructing the following 
spectral expansion, 



v 2 (p) 



c 2 exp 



[11) 



Alternatively, it might be desirable to enforce both the 
thermodynamic stability condition and the "causality" 
conditions, 2 < v 2 < c 2 , by constructing the spectral 
expansion in the following way, 



v 2 (p) 



1 + exp 



(12) 



For the remainder of this paper, spectral representa- 
tions of the equation of state will be based on the fa- 
miliar adiabatic index T(p). The discussion in Sec. IIIII 
shows that T(p) is a reasonably slowly varying function 
for "realistic" neutron-star equations of state, which can 
be represented fairly accurately using expansions having 
only a few terms. Using T(p) for these expansions also 
allows us to make straightforward comparisons with pub- 
lished piecewise-polytrope approximations to the equa- 
tion of state 0. The accuracy and efficiency of these 
expansions in representing "realistic" neutron-star equa- 
tions of state is explored in Sec. IIIII 



B. Enthalpy Based Forms 

The spectral expansions of the standard representa- 
tion of equation of state, e = e(p), should be quite useful 



2 The condition v 2 < c 2 only represents a true causality condition 
if the equation of state e(p) describes both the time dependent 
and the equilibrium properties of the material. In neutron-star 
matter, various strong and weak nuclear interactions determine 
the relative abundances of the various particle species in the equi- 
librium state. Thus v 2 evaluated for the equilibrium equation of 
state only describes the sound propagation speed for low enough 
frequency waves that the material remains continuously in equi- 
librium. The condition v 2 < c 2 may or may not represent a 
causality constraint therefore on sound waves with short enough 
wavelengths to be physically relevant in neutron stars. 



for many applications. For some applications, however, 
the standard representation, e = e(p), is not ideal. For 
example, a useful form of the relativistic stellar struc- 
ture equations Q requires the equation of state to be 
expressed in terms of the relativistic enthalpy, h. For 
applications such as this, e = e(p) must be re-written as 
a pair of equations e = e(h) and p = p{h), where h is 
defined as 

hir)= [wh- (13) 

The needed expressions, e = e(h) and p = p(h), are con- 
structed by inverting h = h(p) from Eq. (|13|) to obtain 
p = p{h), and composing the result with the standard 
equation of state, e = e(p), to obtain e(h) ~ e[p{h)\. 

The transformations needed to construct e = e(h) and 
p = p(h) are difficult to perform numerically in an effi- 
cient and accurate way. Therefore it may be preferable 
to construct a spectral expansion of the equation of state 
based directly on h. This can be done using the methods 
described above for the standard e = e(p) representation. 
To do this a spectral expansion of the adiabatic index, 
considered now as a function of the enthalpy T(h), must 
be defined. The scaled enthalpy variable x = log(h/ho), 
is found to be useful, where ho is the lower bound on the 
enthalpy, ho < h, in the domain where the spectral ex- 
pansions are constructed. The expression for T(x) given 
in Eq. Q then provides a useful expansion for T(h): 



T(h) = exp ■ 



(14) 



k k 

Next, the functions p(h) and e{h) are defined by the sys 
tern of ordinary differential equations, 

dp 

dh = £+P ' 
(te (e+p) 2 

dh ~ pT(h) 



(15) 
(16) 



that follow from the definitions of h, Eq. (|13p . and T, 
Eq. @. The general solution to these equations can be 
reduced to quadrature: 



p(h) = p exp 



e(h) = p{h)- 



h h' 



dh' 



h ~Kh') 

m 



where fl(h) is defined as, 



Po e" 
eo +Po 



r(h>) 



dh'. 



(17) 
(18) 

(19) 



The constants po and eo are defined by po = p{ho) and 
eo = e(h ) respectively. While these quadratures can not 
be done analytically for the spectral expansion defined in 
Eqs. (|14[) . they can be done numerically very efficiently 
and accurately using Gaussian quadrature, as in the stan- 
dard equation of state case. 



III. SPECTRAL FITS TO REALISTIC 
EQUATIONS OF STATE 

The discussion in Sec. HI] shows how any equation of 
state can be represented by spectral expansions of the 
adiabatic index, like the one given in Eq. ((6j) . When these 
expansions are truncated, keeping only a finite number of 
terms, they produce fits to the equation of state, £fit(p) 
that are expected to converge to the exact e(p) as the 
number of terms in the expansion increases. In analogy 
with Fourier series, the rate of convergence for these fits 
should be exponential for smooth equations of state, and 
power law for less than smooth cases. The smoothness of 
an equation of state is determined by the details of the 
micro-physics that controls the properties of the material. 
The convergence rate of an expansion will be reduced 
therefore, from exponential to power law, when phase 
transitions or other non-smooth transitions are present. 
The number of terms required to achieve a certain level of 
accuracy in eat(p), therefore, will depend on the smooth- 
ness and variability of the adiabatic index T(p), and the 
suitability of the chosen spectral basis functions 

The accuracy and practicality of spectral expansions 
for two forms of the equation of state, e = e(p) and 
e = e(/i), are evaluated in this section by construct- 
ing fits to 34 "realistic" neutron-star equations of state. 
These spectral fits arc based on finite spectral expan- 
sions of r(p) and r(/i) respectively. The equations of 
state used for these fits are the same as those used by 
Read, Lackey, Owen and Friedman Q in their study 
of piecewise-polytrope approximations. These realistic 
equations of state are based on a variety of different mod- 
els for the composition of neutron-star matter, and a va- 
riety of different models for the interactions between the 
particle species present in the model material. Descrip- 
tions of these realistic equation of state models, and ref- 
erences to the original publications on each of these equa- 
tions of state are given in Ref. 0, and are not repeated 
here. The individual equations of state are referred to 
here using the abbreviations used in Ref. 0, e.g. PAL6, 
APR1, BGN1H1, etc. The list of these equations of state 
are given in the first column of Table III of Ref. @ , and 
the first columns of Tables HI and [TTI in this paper. 

Approximate equations of state, egt(p), have been con- 
structed for each of the realistic neutron-star equations of 
state listed in Table |TJ These approximations are based 
on Eq. ([7} with T(x) determined by the spectral expan- 
sion in Eq. ([6]). The em(p) constructed in this way de- 
pend on the pressure through the variable x = log(p/po), 
as well as the spectral coefficients 7^: egt = efitO^Tfe)- 
The optimal choice of spectral coefficients, 7/., is made by 
minimizing the differences between emi^i, 7fc) and the ex- 
act ti = e(xi) for a set of pressures, Xi, from the realistic 
neutron-star equation of state tables. These differences 
are measured by constructing the residual: 
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FIG. 1: Residuals Ask are illustrated for several spectral fits 
and for the polynomial fit A p4 for each of the realistic equa- 
tion of state models, which are represented as points along 
the horizontal axis in this figure. These residuals are for stan- 
dard pressure-based representations of the equation of state, 
e — e(p), which are also given in Table U 



where the sum is over all the pressures in the tabulated 
realistic equation of state in the range po < Pi < Pmax- 
These are the pressures that may be present in the cores 
of neutron stars where the equation of state is not well 
known: po is the pressure where the baryon density is 
po = 2 x 10 14 g/cm 3 , and p ma , x is the central pressure 
of the maximum-mass non-rotating neutron-star model 
for the particular equation of state. The constants po, 
log(p max /p ), and the total energy density, e = 
e(po), are given (in cgs units) in the last three columns 
of Table |T] for each of the realistic equations of state. 
This range of pressures coincides with the range used 
by Read, Lackey, Owen and Friedman to construct 
their piecewise-polytrope fits. Using the same range of 
pressures here makes comparison with their work more 
straightforward. 

The spectral coefficients, 7^, that determine the partic- 
ular £fit(a;,7fc) are chosen to minimize the residual A (7*;) 
for each realistic equation of state. The minimization 
process was carried out with an algorithm based on the 
Levenberg-Marquardt method [l||, starting with initial 
estimates, 70 = 1 and jk = for k > 1. Approximate egt 
were constructed in this way for expansions containing 
2, 3, 4, and 5 spectral basis functions. The minimum 
values of the residuals, As2j A53, A 54, and A 55, for 
these cases are listed for each equation of state in Ta- 
ble |TJ The average values of these minimum residuals de- 
crease from about 2.9% for the 2-parameter fits, to about 
0.9% for the 5-parameter fits. For comparison the resid- 
uals Ap4 for the 4-parameter piecewise-polytrope fits of 
Reid, Lackey, Owen, and Friedman Q are also given in 
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TABLE I: Spectral Expansions of the Standard e = e(p) Form of Realistic Neutron-Star Equations of State 



EOS 
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7o 


7i 


72 


73 






Pi 






e< 


/c 2 
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0.0181 
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01 
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2 
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5 


79 
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1 
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X 


10 33 
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6 


73 


API 
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7 


38 
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0091 
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7065 


0.2866 
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34 
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7 


27 
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0204 
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6 


87 
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0.1548 
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0034 
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7 
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0062 
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8079 
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2 
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X 
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7 


25 


WFF3 





0314 





0081 





0081 





0060 





0059 


1 


4126 


-0.1797 


0.0389 


-0.0035 





81 


X 


10 33 


2 


03 


X 


10 14 


7 


32 


BBB2 





1016 





0283 





0238 





0167 





0042 





7390 


0.4555 


-0.1406 


0.0121 


1 


37 


X 


10 33 


2 


05 


X 


10 14 


6 


97 


BPAL12 





0739 





0138 





0086 





0043 





0018 


1 


1081 


-0.3078 


0.0891 


-0.0081 


2 


51 


X 


10 33 


2 


05 


X 


10 14 


6 


08 


ENG 





0527 





0181 





0168 





0138 





0118 





9820 


0.2716 


-0.0862 


0.0075 


1 


33 


X 


10 33 


2 


04 


X 


10 14 


6 


98 


MPA1 





0365 





0223 





0022 





0022 





0019 


1 


0215 


0.1653 


-0.0235 


-0.0004 


1 


51 


X 


10 33 


2 


04 


X 


10 14 


6 


63 


MSI 





0581 





0256 





0052 





0039 





0004 





9189 


0.1432 


0.0122 


-0.0094 


4 


54 


X 


10 33 


2 


04 


X 


10 14 


5 


05 


MS2 





0155 





0074 





0013 





0005 





0001 





9598 


-0.0527 


0.0091 


-0.0035 


4 


10 


X 


10 33 


2 


04 


X 


10 14 


4 


77 


MSlb 





0206 





0179 





0058 





0032 





0004 


1 


2132 


-0.0648 


0.0561 


-0.0111 


3 


18 


X 


10 33 


2 


03 


X 


10 14 


5 


40 


PS 





0568 





0566 





0294 





0284 





0182 


1 


3896 


-0.8472 


0.2636 


-0.0218 


5 


49 


X 


10 33 


2 


05 


X 


10 14 


4 


77 


GS1 





0536 





0762 





0385 





0333 





0265 


1 


8662 


-1.4266 


0.4450 


-0.0389 


3 


22 


X 


10 33 


2 


04 


X 


10 14 


6 


28 


GS2 





0416 





0582 





0427 





0425 





0294 


1 


4580 


-0.7219 


0.1828 


-0.0117 


4 


06 


X 


10 33 


2 


04 


X 


10 14 


4 


99 


BGN1H1 





0435 





0792 





0460 





0439 





0328 


1 


3450 


-0.0996 


-0.0833 


0.0161 


2 


30 


X 


10 33 


2 


05 


X 


10 14 


6 


40 


GNH3 





0092 





0130 





0090 





0081 





0057 


1 


0366 


-0.0044 


-0.0440 


0.0075 


3 


77 


X 


10 33 


2 


06 


X 


10 14 


5 


26 


HI 





0226 





0200 





0117 





0089 





0069 


1 


0653 


0.0362 


-0.1098 


0.0179 


3 


17 


X 


10 33 


2 


04 


X 


10 14 


5 


04 


H2 





0300 





0181 





0133 





0072 





0069 


1 


0743 


0.2250 


-0.2029 


0.0290 


2 


96 


X 


10 33 


2 


04 


X 


10 14 


4 


98 


H3 





0308 





0130 





0109 





0086 





0066 


1 


1340 


0.0925 


-0.1303 


0.0190 


3 


12 


X 


10 33 


2 


04 


X 


10 14 


4 


91 


H4 





0098 





0098 





0097 





0069 





0063 


1 


0526 


0.1695 


-0.1200 


0.0150 


3 


16 


X 


10 33 


2 


04 


X 


10 14 


5 


21 


H5 





0214 





0150 





0126 





0054 





0054 


1 


0106 


0.2765 


-0.2011 


0.0270 


2 


97 


X 


10 33 


2 


04 


X 


10 14 


5 


09 


H6 





0185 





0137 





0133 





0130 





0100 


1 


0650 


-0.0196 


-0.0474 


0.0077 


3 


35 


X 


10 33 


2 


04 


X 


10 14 


4 


80 


H7 





0139 





0132 





0103 





0059 





0056 





9582 


0.1619 


-0.1294 


0.0177 


3 


19 


X 


10 33 


2 


04 


X 


10 14 


5 


22 


POL2 





0227 





0252 





0121 





0090 





0075 


1 


0410 


0.0173 


-0.0904 


0.0150 


2 


90 


X 


10 33 


2 


04 


X 


10 14 


5 


44 


ALF1 





0947 





0669 





0453 





0369 





0305 


1 


0143 


-0.3102 


0.1809 


-0.0248 


1 


50 


X 


10 33 


2 


05 


X 


10 14 


6 


21 


ALF2 





0655 





0629 





0450 





0256 





0230 





4613 


1.5237 


-0.5817 


0.0571 


1 


50 


X 


10 33 


2 


05 


X 


10 14 


5 


79 


ALF3 





0371 





0355 





0139 





0132 





0131 





8536 


0.2405 


-0.0743 


0.0041 


1 


50 


X 


10 33 


2 


05 


X 


10 14 


6 


14 


ALF4 





0453 





0652 





0166 





0089 





0088 





8806 


0.0656 


0.0765 


-0.0177 


1 


50 


X 


10 33 


2 


05 


X 


10 14 


5 


97 


Average 


0.0383 


0.0287 


0.0149 


0.0114 


0.0085 







Table Q] for each equation of state. 3 The values of the 
residuals, A52, As3, Ag4, A55, and Ap 4 , are also shown 



3 The residuals reported in Table III of Ref. [i| differ from the 
Ap^ residual listed in Table J] in two ways. The first difference 
is the residuals reported in Ref. Q| are evaluated using base- 
10 logarithms, rather than the natural logarithms used here. 
This difference makes the Ref. Q| residuals smaller by the factor 
log 10 e ?s 2.3. The second difference is that Ap4 reported here 
measures the accuracy of the piecewisc-polytrope fits for e(p), 
while the residuals reported in Ref. Q measure the accuracy of 
those fits for p(p) where p is the baryon density of the material. 



graphically in Fig. [T] Points along the horizontal axis in 
Fig. [T] represent the different realistic equations of state 
in the order listed in Table Q] These results show that 
the spectral fits are convergent, and do a fairly good job 
of approximating this collection of realistic equations of 
state. The 2-parameter spectral fits have smaller residu- 
als than the 4-parameter piecewise-polytrope fits for most 
of these equations of state. Also listed in Table fl] are the 
optimal values of the spectral coefficients, 70, 71, 72, and 
73, for the 4-parameter spectral approximation to each 
equation of state. These values, together with the tabu- 
lated constants eo and po can be used to re-construct the 
complete 4-paramcter spectral fits for e(p), using Eqs. © 
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2.4 



2.3 



2.2 



2.1 



PAL6 eos data 
4-param polytropic fit 

4-param spectral fit 




1 2 3 4 5 6 

x = log(p/p ) 

FIG. 2: Adiabatic index as a function of pressure, r(x), for 
various fits to the PAL6 equation of state. 







i 1 i 1 i 1 i 1 i - ■ — r 

• BGN1H1 eos data 
4-param polytropic fit 

— 4-param spectral fit 



2 3 4 5 
x = log(p /p Q ) 



FIG. 4: Adiabatic index as a function of pressure, T(a;), for 
various fits to the BGN1H1 equation of state. 



r 



• APRl eos data 

— 4-param polytropic fit 

- 4-param spectral fit 




./ 

7 



2 4 6 8 

x = log(p/p ) 

FIG. 3: Adiabatic index as a function of pressure, T(a;), for 
various fits to the APRl equation of state. 



0.02 



w 



co 0.01 - 



&J0 

o 



PAL6 



4-param polytropic fit 

2- param spectral fit 

3- param spectral fit 

4- param spectral fit 

5- param spectral fit 




1 2 3 4 5 
x = log(p/p Q ) 

FIG. 5: Ratios between various fits, eat(x), and the exact 
PAL6 equation of state, e(x). Note that log(eflt/e) « (eat — 
e)/e measures the fractional error of the fit. 



and 0. 

Three equations of state have been chosen from the 
complete set to illustrate in more detail the accuracy of 
the fits. These three cases are equation of state PAL6 
having the highest accuracy spectral fits, APRl having 
spectral fits with average accuracy, and BGN1H1 having 
the lowest accuracy spectral fits. Figures [20 show the 
adiabatic index T(x) computed directly from the tab- 
ulated equations of state, the 4-parameter spectral fit 
to the adiabatic index, and the 4-parameter piecewise- 
poly trope fit for each of these equations of state. From 



these figures it is clear that the smoother equations of 
state, like PAL6, have the highest accuracy fits, while 
the equations of state with a sharp phase transition, like 
BGN1H1, have the lowest accuracy fits. Figures [5][7] illus- 
trate the errors in the various fits to the equation of state 
e(x) for these three example equations of state. These er- 
rors are illustrated as graphs of log(efit/e) ~ {e&t — e)/e, as 
functions of the pressure variable x = \og(p/p ). These 
figures show that the spectral fits are much more ac- 
curate than the 4-parameter piecewise-polytrope fit for 
the smooth and average equations of state, PAL6 and 




FIG. 6: Ratios between various fits, e^ t {x), and the exact 
APR1 equation of state, e(x). Note that log(eflt/e) ~ (em — 
e)/e measures the fractional error of the fit. 



0.3 - 



BGN1H1 



1 - 1 — i — 1 — i — 1 — r 

4-param polytropic fit 

2- param spectral fit 

3- param spectral fit 

4- param spectral fit 

5- param spectral fit 




1 2 3 4 5 6 
x = log(p/p ) 



FIG. 7: Ratios between various fits, ea t (x), and the exact 
BGN1H1 equation of state, e(x). Note that log(efitA) ~ (est — 
e)/e measures the fractional error of the fit. 



APR1 . The spectral fits have about the same accuracy as 
the piecewise-polytrope fit for equation of state BGN1H1 
which has a strong phase transition. Figures [5][7] also il- 
lustrate in a visual way the convergence of the spectral 
fits as the number of basis functions is increased. 

The enthalpy based representation of the equation of 
state, e = e(h) and p = p(h), is more useful for certain 
purposes. It is helpful to understand, therefore, whether 
the spectral fits for this representation, Eqs. (fT4")) -(fTi? |) . 



FIG. 8: Residuals Ask are illustrated for several spectral fits 
and for the polynomial fit A p4 for each of the realistic equa- 
tion of state models, which are represented as points along the 
horizontal axis in this figure. These residuals are for enthalpy- 
based representations of the equation of state, e = e(h), which 
are also given in Table [Til 



are as accurate and effective as those for the standard 
e = e(p) representation. So approximate equations of 
state, efit(^), have been constructed for the same set of 
34 realistic neutron star equations of state described in 
Ref. 0. As before these fits were made by adjusting the 
values of the spectral coefficients jk defined in Eq. (|I4|) . 
to minimize the residual A (7*;) defined in (f2"U)) . The only 
differences between this and the standard case are: The 
variable x = log(/i//io) is chosen here to be an enthalpy 
variable, and the functions e(h,"/k) and p(h,jk) are de- 
termined here with Eqs. ((TT1) — (JTHJ) . Table [TT] lists each 
equation of state, along with the residuals, Ag2, Ag3, 
Ag4, and Ag4, for the spectral fits having 2, 3, 4 and 5 
non-zero spectral coefficients. The values of the residuals, 
A52, As3, As4, A55, and Ap4, are also show graphically 
in Fig. [5] for these enthalpy based fits. Table HTl also lists 
the coefficients, 70, 71, 72, and 73 for the 4-parameter 
spectral fit. Finally, Table [TT1 contains information about 
the enthalpy variables, ho and x max = fog(/i m ax/^o)j 
for each equation of state. The quantity ho is the en- 
thalpy for which eo = e(ft.o) and po = p(ho) whose values 
are listed in Table [H and h max is the value for which 
Pmax = p{h m a,x)- These results show that the spectral 
expansions of the enthalpy based representation of the 
equation of state are very comparable to the standard 
expansions. 
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TABLE II: Spectral Expansions of the e = e(h) Form of Realistic Neutron-Star Equations of State 



EOS 




A S2 




A S 3 


A S4 


A S5 




7o 


7i 


72 


73 




h 


3^ max 


PAL6 





0032 





0016 





0005 





0002 





8608 


-0.1509 


0.0909 


-0.0192 





0405 


2 


62 


SLy 





0089 





0035 





0017 





0006 


1 


0077 


0.2084 


-0.1266 


0.0203 





0314 


3 


22 


API 





0708 





0331 





0186 





0081 





5205 


1.3288 


-0.7804 


0.1316 





0276 


3 


31 


AP2 





0307 





0139 





0063 





0034 





7557 


0.6423 


-0.3586 


0.0604 





0287 


3 


34 


AP3 





0250 





0054 





0024 





0015 





9520 


0.4691 


-0.2109 


0.0273 





0298 


3 


40 


AP4 





0249 





0030 





0014 





0014 





8824 


0.4064 


-0.1405 


0.0135 





0297 


3 


41 


FPS 





0044 





0041 





0037 





0018 


1 


1462 


-0.1088 


0.0432 


-0.0106 





0279 


3 


23 


WFF1 





0521 





0105 





0068 





0064 





7115 


0.8263 


-0.3258 


0.0372 





0273 


3 


53 


WFF2 





0267 





0118 





0070 





0037 





8527 


0.6162 


-0.3135 


0.0512 





0285 


3 


48 


WFF3 





0109 





0101 





0058 





0054 


1 


3660 


-0.4236 


0.2293 


-0.0469 





0253 


3 


30 


BBB2 





0318 





0277 





0210 





0089 





8447 


0.8101 


-0.5808 


0.1120 





0297 


3 


24 


BPAL12 





0168 





0097 





0048 





0022 


1 


0847 


-0.6537 


0.4289 


-0.0874 





0373 


2 


73 


ENG 





0195 





0179 





0142 





0109 


1 


0426 


0.4716 


-0.3353 


0.0632 





0286 


3 


39 


MPA1 





0251 





0032 





0030 





0022 


1 


0523 


0.3546 


-0.1356 


0.0074 





0297 


3 


33 


MSI 





0277 





0055 





0035 





0003 





9340 


0.2231 


0.0718 


-0.0642 





0485 


2 


73 


MS2 





0096 





0021 





0003 





0001 





9680 


-0.1326 


0.0786 


-0.0360 





0367 


2 


56 


MSlb 





0192 





0062 





0029 





0003 


1 


2148 


-0.1255 


0.1959 


-0.0700 





0403 


2 


92 


PS 





0624 





0298 





0298 





0166 


1 


3064 


-1.2479 


0.5534 


-0.0186 





0571 


2 


26 


GS1 





0889 





0383 





0377 





0249 


1 


7649 


-1.9609 


0.8871 


-0.0942 





0265 


3 


22 


GS2 





0603 





0437 





0433 





0314 


1 


2562 


-0.4936 


-0.1125 


0.1051 





0362 


2 


61 


BGN1H1 





0868 





0495 





0439 





0403 


1 


1915 


0.2463 


-0.7222 


0.2300 





0356 


2 


97 


GNH3 





0135 





0093 





0081 





0057 


1 


0236 


0.0297 


-0.2197 


0.0742 





0482 


2 


53 


HI 





0200 





0122 





0084 





0071 


1 


0365 


0.1735 


-0.5519 


0.1790 





0406 


2 


39 


H2 





0173 





0136 





0067 





0067 


1 


0666 


0.4639 


-0.8300 


0.2393 





0391 


2 


41 


H3 





0124 





0109 





0082 





0068 


1 


1204 


0.2084 


-0.5285 


0.1519 





0402 


2 


42 


H4 





0103 





0100 





0066 





0066 


1 


0070 


0.4040 


-0.3709 


0.0696 





0223 


3 


18 


H5 





0137 





0134 





0059 





0049 





9614 


0.6257 


-0.6459 


0.1371 





0233 


2 


96 


H6 





0137 





0135 





0129 





0113 


1 


0182 


0.1483 


-0.2638 


0.0628 





0263 


2 


79 


H7 





0123 





0112 





0057 





0057 





9147 


0.4331 


-0.4779 


0.1064 





0254 


2 


94 


FOL2 





0265 





0130 





0086 





0077 


1 


0127 


0.1289 


-0.4822 


0.1616 





0389 


2 


52 


ALF1 





0731 





0475 





0400 





0283 





9349 


-0.3810 


0.6697 


-0.2305 





0305 


2 


75 


ALF2 





0692 





0490 





0276 





0178 





7100 


2.6577 


-2.2438 


0.4712 





0305 


2 


83 


ALF3 





0386 





0149 





0140 





0134 





8987 


0.4588 


-0.3533 


0.0489 





0305 


2 


73 


ALF4 





0708 





0135 





0099 





0094 





8705 


0.4469 


0.0284 


-0.0914 





0305 


2 


83 


Average 


0.0323 


0.0166 


0.0124 


0.0089 







IV. DISCUSSION 

The spectral fits constructed here were designed to ex- 
plore how accurately the real neutron-star equation of 
state might be determined once the first few accurate 
neutron-star mass-radius measurements become avail- 
able. These fits could be improved for equations of state 
with phase transitions by using separate spectral fits 
above and below the phase-transition pressure. These 
piecewise spectral fits could eliminate Gibbs phenomena 
errors, but would require significantly more (roughly dou- 
ble the number of) parameters. Measuring this larger 
number of parameters from neutron-star observations 



would therefore require a much larger number of accu- 
rate mass-radius determinations. A systematic study of 
the relative accuracy of single versus piecewise spectral 
expansions (with the same total number of free param- 
eters) will therefore be delayed until more neutron-star 
measurements become available. 

The fits constructed here show that spectral represen- 
tations of the various realistic neutron-star equations of 
state are remarkably accurate, even when the number of 
basis functions used in the spectral expansion is rather 
small. Such expansions provide an attractive alternative 
to the piecewise-polytrope approximations for a variety of 
current "realistic" models of the equation of state. If the 
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real neutron-star equation of state is relatively smooth, 
then these spectral expansions will provide an extremely 
efficient way to represent it. And remarkably, even if 
the neutron-star equation of state has a phase transition, 
these spectral fits do about as well as the piecewise poly- 
trope fits with the same number of free parameters. 
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